Inhomogeneity, Dynamical Symmetry, and Complexity in High- Temperature Superconductors: 
Reconciling a Universal Phase Diagram with Rich Local Disorder 
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A model for high- temperature superconductors incorporating antiferromagnetism, rf-wave superconductivity, 
and no double lattice-site occupancy can give energy surfaces exquisitely balanced between antiferromagnetic 
and superconducting order for specific ranges of doping and temperature. The resulting properties can reconcile 
a universal cuprate phase diagram with rich inhomogeneity, relate that inhomogeneity to pseudogaps, give a fun- 
damental rationale for giant proximity effects and other emergent behavior, and provide an objective framework 
for separating essential from peripheral in the superconducting mechanism. 
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High-temperature superconductivity was discovered more 
than two decades ago [1], but its interpretation remains con- 
troversial |2]. We believe that this lack of theoretical consen- 
sus may be attributed to two fundamental, and related, issues: 
(1) Most models emphasize a limited aspect of the (complex) 
problem; as a result, there are few solvable models that cap- 
ture a sufficient range of essential physics. (2) The complex 
behavior of these compounds obscures the superconducting 
mechanism because, in the absence of solvable models incor- 
porating a wide enough range of physics, it is difficult to sep- 
arate essential features from secondary ones using only data. 

For example, high-temperature superconductors exhibit a 
variety of spatial inhomogeneities such as stripes or checker- 
boards, particularly in the hole-underdoped region and near 
magnetic vortex cores lH i, 11 H H . The relationship of this 
inhomogeneity to the unusual properties of these systems is 
unsettled. Does it lead to superconductivity (SC), does it op- 
pose SC, or is it a sideshow? Particularly vexing is that this 
rich variety of inhomogeneity would at first glance seem to 
be inconsistent with broad evidence for a cuprate phase dia- 
gram that is rather universal, particularly for hole-doped com- 
pounds. How does one reconcile a seemingly universal phase 
diagram with a bewildering aiTay of inhomogeneity for indi- 
vidual compounds? 

Dopant atoms play a dual role in high-temperature super- 
conductors: they support SC globally by enhancing charge 
carrier density, but may suppress SC locally through atomic- 
scale disorder. McElroy et al |8] found strong disorder in 
atomically-resolved scanning tunneling microscope images of 
the SC gap for Bi-2212. They concluded that this disorder de- 
rives primarily from dopant impurities and that charge varia- 
tion between nanoregions is small, implying that inhomogene- 
ity is tied to impurities need not couple strongly to charge. 

We show here that such inhomogeneities are a generic con- 
sequence of perturbations on the antiferromagnetic (AF) and 
SC correlations, largely independent of specifics and not nec- 
essarily coupled to charge variation. Further, we show that 
these properties are consistent with a global cuprate phase di- 
agram, are tied intimately to the nature of pseudogap states, 
and imply a linkage among pseudogaps, inhomogeneity, and 
emergent behavior Thus we provide a testable hypothesis for 



separating primary from derivative features in the rich high- 
temperature superconductor data set. 

The SU(4) model |9, 10, 1 1, 12, iKIl is a fermion many- 
body theory that incorporates AF and SC order on an equal 
footing, conserves spin and charge, and implies no double oc- 
cupancy on the lattice (Mott insulator ground state at half fill- 
ing). The effective SU(4) Hamiltonian can be expressed as 

H = Ho- Go [(1 - a)D'D + aQ ■ Q] +8'S-S, (1) 

where Hq, Go and g' are constants, creates cZ-wave sin- 
glet pairs, Q is the staggered magnetization, S is spin, a = 
a{x) = x(x)/(x(jc) -HGo(jc)), and Go = x{x)+Go{x), where 
Go{x) and are effective SC and AF coupling strengths. 
Doping is characterized by a parameter .it = 1 —n/Q. for an n- 
electron system, with Q. the maximum number of doped holes 
(or doped electrons for electron-doped compounds) that can 
form coherent pairs, taking the normal state at half filling as 
the vacuum. Generally, P ~ 0.25ji:, where P is the standard 
hole-doping parameter, normalized to the number of copper 
sites [lOfl. 

Three dynamical symmetry limits have exact solutions i^. 
The SO(4) limit (a = 1) is a collective AF state, the SU(2) 
limit (a = 0) is a collective SC state, and the SO(5) limit 
(a = 1 /2) is a critical dynamical symmetry interpolating be- 
tween the SC and AF limits |j9t]. For other values of cr an 
approximate solution can be obtained using generalized co- 
herent states. We characterize AF in these solutions by the 
staggered magnetization Q = {Q^) , or hy fi defined through 



e = 2i2j3(«/(2£2)-j3'^) 



2n1/2 



(2) 



and singlet d-wscve pairing through A = (D^D) . 

Total energy surfaces are obtained from the ejcpectation 
value of Eq. ([T]) in coherent state approximation |9]. To re- 
late them to data we use the variation of a with doping, which 
was determined by fitting to cuprate data and given in Fig. 1 of 
Ref. 1 1 l]a. Figure[T]shows SU(4) energy surfaces as functions 
of AF order j3, SC order parameter A, and doping x. These 
energy surfaces exhibit two fundamental instabilities that may 
play a large role in the properties of cuprate superconductors. 
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FIG. 1: (Color online) Total energy vs. (a) AF order /3 and (b) SC or- 
der A; curves labeled by hole-doping x ~ 4P, where P is the number 
of holes per copper site. Energy is in units of 1^, with x the AF 
coupling strength. The black dashed line indicates the critical doping 
X = Xq (see Ref. 1 11]); color coding indicates energy surfaces favor- 
ing SC (red), AF (blue), and AF -I- SC (green). The doping-dependent 
ranges of the energy contours in the order parameters reflect the finite 
valence space (single band) of the model. 



The first is an instability against condensing pairs when the 
system is doped away from half-filling infinitesimally. This 
instability is similar to the Cooper instability for a normal su- 
perconductor, but generalized to a doped Mott insulator, and 
accounts for the remarkably rapid development of supercon- 
ductivity with hole doping in the cuprates [12J. In this paper 
we concentrate on the properties and implications of a second 
fundamental instability that occurs in the underdoped region. 

The energy surfaces at constant doping fall into three gen- 
eral classes: AFh-SC (e.g., x = 0.1), SC (e.g., x = 0.9), and 
critical (e.g., x ~ 0.6, which marks a quantum phase transi- 
tion). Curves in the AFh-SC class have minima at finite and 
large j3o, and small but finite Ao, where the subscript zero de- 
notes the value of the order parameter at the minimum of the 
energy surface. Curves in the class SC are characterized by 
j3o = and finite Aq. Of most interest here are surfaces near 
critical in Fig.IT] which correspond to broken SU(4) D SO(5) 
dynamical symmetry |@] and are flat over large regions of pa- 
rameter space. This implies that there are many states lying 
near the ground state with very different values for j3 and A. 
Thus the surface is critically balanced between AF and SC or- 
der, and small perturbations can drive it from one to the other 
This defines a critical dynamical symmetry of the SU(4) alge- 
bra J^; we shall term this situation dynamical criticality. 

The extreme sensitivity of critical surfaces to perturbations 
is illustrated in Fig. |2l Each set of curves is associated with 
a fixed value of doping x = 0.6 (equivalently, P = 0. 15), with 
the solid line corresponding io a — 0.6, the dashed line to a 
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FIG. 2: (Color online) Energy surfaces as a function of the AF order 
parameter jS for fixed doping P = x/4 = 0.15. The dashed line cor- 
responds to (7 = cTq = 0.6, the upper curve to a 10% increase in ct, 
and the lower curve to a 10% reduction in CT. 



10% increase in a (AF perturbation), and the dotted line to a 
10% reduction in a (SC perturbation). The effect on the en- 
ergy surface versus A (not shown) is significant but less dra- 
matic: Aq is shifted, but remains finite in all three cases. We 
see that this small fluctuation in a can alter the energy surface 
between AFh-SC (finite )3o and Aq) and SC (/3o = and finite 
Ao). This sensitivity is specific to the critical (broken SU(4) D 
SO(5)) dynamical symmetry. The AF region near x = and 
the (i-wave superconducting region at larger hole doping (see 
Fig.[T|) are very stable against such perturbations. 

The AF instability displayed graphically in Fig.|2]may also 
be understood analytically. From the T = solution for Q 
given by Eqs. (24a) and (14) of Ref. | ll]a, we find 



dx 



1 



' 2x 



^[{Xq-x){Xq'~x)Yl^ 



(3) 



and a small change in doping will cause a large change in anti- 
ferromagnetic correlations nearx Xq. This is a consequence 
of SU(4) symmetry, which requires that Q vanish for x > x^ 
and be finite for < x < x^. 

It is not hard to conjecture mechanisms altering the ratio of 
AF to SC coupling locally. For example, Ref. [8] found that 
nanoscale disorder is tied to influence of dopant impurities. 
Nunner et al [15] (see also [16]) compared these results with 
Boboliubov-de Gennes calculations and proposed that out-of- 
plane dopant atoms can modulate pairing on a scale compara- 
ble to the lattice spacing, through lattice-distortion modifica- 
tion of electron-phonon coupling or superexchange. 

States associated with critical dynamical symmetry may be 
termed chameleon states: their variational energy surfaces 
are flat over large regions of parameter space and their in- 
trinsic collective properties may be changed qualitatively by 
a small perturbing background that alters the AF-SC compe- 
tition. Figure [T] suggests that underdoped cuprates have near- 
critical energy surfaces. Thus, chameleon states are central to 
the discussion of inhomogeneity and to the general issue of 
understanding pseudogap states in underdoped cuprates. 

Figure [3] is constructed from the expectation value of ([T]| 
in coherent state approximation, assuming a 1-D spatial per- 
turbation, sin(2;rL), with a = 0.6. It illustrates schematically 
how a small (10%) periodic fluctuation in the AF and SC cou- 
pling for a critically symmetric underdoped compound can 
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FIG. 3: (Color online) A small, spatially-periodic variation in the 
AF and SC coupling can produce inhomogeneity without significant 
charge variation. The charge is M and a subscript m denotes maxi- 
mum values. 



lead to inhomogeneity. In this example, 1 -dimensional spatial 
variations of the coupling ratio a give fluctuations in order 
parameters leading to stripes in which AFh-SC (c7 > 0.6) and 
SC (a < 0.6) are favored alternately. Also shown are the re- 
sponses of AF fluctuations dfi /dL to this variation in a. (We 
do not intend this as a realistic model of a stripe phase, but as 
a cartoon indicating how such a model could be built.) 

Figure[3]indicates that a small spatial variations in the cou- 
pling ratio (7 can produce regions having large AF and weaker 
SC order, interspersed with regions having significant SC cor- 
relations but no AF order. Although not plotted, the AFh-SC 
regions also exhibit small triplet pair densities, which vanish 
in the SC regions. The primary fluctuation between stripes is 
in the AF order, which can jump between zero and the maxi- 
mum allowed by the SU(4) model between adjacent stripes; 
the pairing order changes are much smaller and pairing is 
finite in both the SC and SCh-AF regions. The variation of 
dp /dL indicates that appreciable softness in AF and SC may 
occur on the boundaries between regions. 

The expectation value of the charge M is not a function 
of the coupling ratio a in the SU(4) coherent state, so the 
critical energy surface fluctuations responsible for alternating 
AFh-SC and SC stripes in Fig. [3] cause no charge variation 
(AM — 0). Data indicate that the relative charge variation for 
Bi-2212 surface nanoscale patches is less than 10%, implying 
heterogeneity not strongly coupled to charge |8]. This view 
is supported by analyses of heat capacity and NMR data on 
Bi-2212 and YBCO that find a universal phase behavior for 
cuprates, with little static charge modulation |.17. 18]. 

Of course, the variation in a could itself be due to a charge 
modulation. From the preceding discussion, we may expect 
that if a charge modulation occurs in either the AF region near 
half filling or the SC region at larger doping, its effect will 
be small because the energy surfaces are not critical there. 
However, if a charge modulation occurs in the underdoped 
region where the energy surfaces are near critical its effect 
could be amplified by dynamical criticality even if a is not 



altered significantly, since this is equivalent to a doping mod- 
ulation. Thus, we suggest a mechanism for producting strong 
inhomogeneity without necessarily invoking charge fluctua- 
tions as the cause, but that can in particular cases for under- 
doped compounds be produced by a charge modulation. Such 
a mechanism could make understandable strong inhomogene- 
ity in the face of a universal phase diagram, and resolve com- 
peting experimental claims regarding the role of charge varia- 
tion in producing inhomogeneity. 

As an aside, inhomogeneity without charge modulation has 
an analog in astronomy. We might assume, erroneously, that 
density in the bright arms of a spiral galaxy (a form of "stripe 
order") is much larger than that between the arms. Actu- 
ally, spiral arms are prominent because the mass (gravitational 
charge) there is more visible: galactic dynamics generates 
many hot, luminous stars in the spiral arms. In analogy, the 
structure in Fig. |3]is not due to charge variation but rather to 
modulation of the charge "visibility" by (quantum) dynamics. 

The minimal patch size that can support SU(4) coherent 
states is crucial to the present argument. Experience (e.g., 
Ref. fW]) suggests that dynamical symmetry can be realized 
in fermion valence spaces having as few as several particles. 
This would be consistent with inhomogeneities on scales com- 
parable to atomic dimensions, as required by data (§]. 

Inhomogeneity caused by electronic self-organization is of- 
ten contrasted with that caused by dopant impurities. Our dis- 
cussion implicates both as sources of nanoscale structure. The 
proximate cause may be impurities that perturb the SU(4) en- 
ergy surface but the criticality of that surface, which greatly 
amplifies the influence of impurities, results from the self- 
organizing, doped Mott insulator encoded in the SU(4) alge- 
bra. Note a recent analysis |16] suggesting, from a different 
perspective, that intrinsic amplification of impurity effects is 
required to explain nanoscale structure in Bi-2212. 

Critical dynamics may also produce inhomogeneities near 
magnetic vortices and magnetic impurities. A magnetic field 
should suppress SC relative to AF, so distance from a vortex d 
may be expected to alter the average value of a, just as chang- 
ing the doping P would. Therefore, we may expect a region 
near magnetic vortices or impurities where the symmetry is 
critical and exhibits sensitivity to perturbations similar to that 
exhibited in Fig.|2] with d modulating a rather than P. 

Giant proximity effects are observed in the cuprates where 
non-SC copper oxide material sandwiched between supercon- 
ducting material can carry a supercurrent, even for a thickness 
much larger than the coherence length [20]. Phenomenology 
indicates that pre-existing nanoscale SC patches can precipi- 
tate such effects ll2lll . We find similar possibiUties on funda- 
mental grounds, but also suggest that the inhomogeneity need 
not pre-exist. Dynamical criticality renders even a homoge- 
neous pseudogap phase unstable against large fluctuations in 
the AF and SC order. Thus, proximity of superconducting ma- 
terial to pseudogap material coupled with perturbations from 
background impurity fields can trigger dynamical nucleation 
of nanoscale structure and giant proximity effects, even if no 
static inhomogeneity exists beforehand. 

In experiments with one unit cell thickness La2Cu04 
AF barrier layers between superconducting Lai.85Sro.i5Cu04 
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samples, Bozovic et al ll22ll found that the two phases did not 
mix, with the barrier layer completely blocking a supercur- 
rent. These results were interpreted to rule out many models 
of high-temperature superconductors, in particular models in 
which SC and AF phases are nearly degenerate like the SO(5) 
model [23]. The absence of a proximity effect between the 
AF and SC phases (but presence of a strong proximity effect 
between pseudogap and SC material) is plausibly consistent 
with the SU(4) model: the AF phase is not rotated directly 
into the SC phase but rather evolves with increased doping 
into a mixed SC and AF phase, which then is transformed into 
a pure SC state at a critical doping point near optimal doping, 
where the AF correlations vanish identically [11]. 

The spontaneous appearance of properties that do not pre- 
exist in a system's elementary components is termed emer- 
gence; systems with emergent behavior are said to exhibit 
complexity (see Dagotto ll24ll for a review). Complexity can 
occur when the choice between potential ground states is sen- 
sitive to even weak external perturbations. The amplification 
effect implied by SU(4) dynamical criticality can facilitate 
emergent behavior and complexity. The giant proximity effect 
and the perturbatively-induced structure of Fig. |3]are exam- 
ples. More generally, critical dynamical symmetry may rep- 
resent a fundamental organizing principle for complexity in 
strongly-correlated fermi systems. For example, critical dy- 
namical symmetries are known in nuclear physics [25]. 

The SU(4) coherent state method that yields the variational 
energy surfaces discussed here admits quasiparticle solutions 
that generalize the BCS equations, giving a rich, highly uni- 
versal phase diagram in agreement with much available data 
ifllll . and a quantitative description of fermi arcs ifisll . There- 
fore, the propensity of the pseudogap state to a broad vari- 
ety of induced inhomogeneity, and a quantitative model of the 
cuprate phase diagram (including pseudogaps) that exhibits 
highly universal character, both follow directly from a model 
that implements AF and c/-wave SC competition in a doped 
Mott insulator This natural coexistence of a universal phase 



diagram with a rich susceptibility to disorder could reconcile 
many seemingly disparate observations in the cuprate super- 
conductors. For example, since the SU(4) pseudogap state has 
both pairing fluctuations and critical dynamical symmetry, it 
could support Nernst vortex states as perturbations on a ho- 
mogeneous phase having incipient nanoscale heterogeneity. 

In summary, competing antiferromagnetism and d-wave 
pairing, constrained by no double site occupancy, can lead to 
energy surfaces in hole-underdoped cuprates that are critically 
balanced between antiferromagnetism and superconductivity. 
These surfaces can be flipped between dominance of one or 
the other by small changes in the ratio of the antiferromag- 
netic to pairing strength. Therefore weak perturbations in the 
underdoped region, and near vortex cores or magnetic impuri- 
ties, can produce amplified inhomogeneity having the spatial 
dependence of the perturbation but the intrinsic character of an 
SU(4) symmetry. (The symmetry defines the possible states; 
the perturbation selects among them.) Our results show that 
such effects could, but need not, imply spatial modulation of 
charge. More generally, we have suggested that critical dy- 
namical symmetry may be a fundamental principle of emer- 
gent behavior in correlated fermion systems. 

The generality of our solution implies that any realistic the- 
ory describing superconductivity competing with antiferro- 
magnetism should contain features similar to those discussed 
here. In particular, similar phenomena may be possible in the 
new iron-based high-temperature superconductors |26]. Then 
the existence of complex inhomogeneities for compounds 
having (paradoxically) universal phase diagrams suggests that 
(1) properties of superconductors in the pseudogap state, and 
near magnetic vortices and impurities, are largely determined 
by critical dynamical symmetry, and (2) inhomogeneity is a 
strong diagnostic for the mechanism of high-temperature su- 
perconductivity but it is only consequence, not cause. 

We thank Takeshi Egami, Elbio Dagotto, Adriana Moreo, 
Pengcheng Dai, and Boris Fine for discussions. 
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